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NOMENCLATURE 


A 

mean  plan  area  of  ice  floes 

P* 

shaft  power 

B 

ship  beam 

Q 

propeller  torque 

CD 
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R 

total  resistance  in  ice 

d 

propeller  diameter 

Si 

Pi/P 

D 

ship  draft 

Sq 

torque  standard  deviation 

E 

ice  modulus 
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F  n 
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fi 
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width  of  ice  floe 
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gravity 
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ship  block  coefficient 

hi 

ice  thickness 

y 

water  specific  weight  ( pg ) 

J 

advance  coefficient  ( Vlnd ) 

A 

ship  displacement 

In 

Atkins  “ice”  number 

X 

geometric  scale 

k\ ,  ho 

coefficients 

K 

velocity  scale 

kq 

torque  coefficient 

% 

propulsion  efficiency 

Kr 

thrust  coefficient 

M 

water  dynamic  viscosity 

Klc 

ice  fracture  roughness 

u 

ice  Poisson’s  ratio 

lc 

ice  characteristic  length 

P 

water  density 

L 

ship  length 

Pi 

ice  density 

Lc 

length  of  ice  floe 

<*c 

ice  compressive  strength 

n 

propeller  rotational  speed 

Of 

ice  flexural  strength 

P 

porosity 

o8 

ice  shear  strength 
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Ship  Model  Testing  in  Level  Ice:  An  Overview 

JEAN-CLAUDE  P.  TATINCLAUX 


INTRODUCTION 

Model  testing  in  ice  of  ships  and  other 
structures  is  a  relatively  recent  branch  of 
marine  engineering,  with  the  first  refriger¬ 
ated  towing  tank  being  built  in  1955  at  the 
Arctic  and  Antarctic  Research  Institute  of  the 
Soviet  Union  in  Leningrad.  However,  with 
the  rapid  development  of  the  arctic  and  sub¬ 
arctic  regions  of  the  globe,  spurred  by  the  dis¬ 
covery  of  large  mineral  and  petroleum  depos 
its,  the  fields  of  arctic  engineering,  in 
general,  and  ice  engineering,  in  particular, 
have  rapidly  grown  in  the  last  decade  or  so. 
For  economic  and  geopolitical  reasons,  the 
number  of  ships  transiting  ice-infested  wa¬ 
ters  has  greatly  increased.  Correspondingly, 
better  and  new  ship  designs  specifically 
adapted  to  the  requirements  of  a  particularly 
harsh  environment  led  to  the  development 
and  construction  of  new  facilities  devoted  to 
ship  model  testing  in  ice.  In  the  two  decades 


from  1955  to  1975,  only  three  ice  towing  tanks 
were  built,  one  in  Finland,  one  in  West  Ger¬ 
many,  the  other  in  the  United  States.  In  the 
following  decade  (1975-1985),  ten  such  tanks 
were  built  in  Japan,  Northern  Europe  (in¬ 
cluding  the  USSR),  and  North  America. 
Table  1  recapitulates  the  ice  tanks  currently 
operating,  with  their  major  characteristics. 

All  types  of  ships  with  varying  missions 
and,  thus,  requirements  are  encountered  in 
the  Arctic  and  Antarctic,  the  Northern  Atlan¬ 
tic,  the  Baltic  sea,  and  as  far  south  as  the 
Great  Lakes  and  St.  Lawrence  seaway;  they 
range  from  tankers  (oil  and  liquified  natu¬ 
ral  gas),  to  cargo  and  container  ships,  supply 
ships,  research  vessels,  and  national  securi¬ 
ty  or  Coast  Guard  vessels. 

As  in  conventional  towing  tank  tests,  the 
main  purposes  of  ship-model  resistance  and 
propulsion  tests  in  ice  tanks  are:  1)  to  predict 
ship  performance  in  ice  for  a  given  design,  2) 
to  suggest  modifications  to  a  proposed  design, 


Table  1.  Existing  ice  testing  facilities.  The  basin  length  given  may  include 
the  length  of  the  trim  tank. 


Main  dimensions 

L  X  B  x  D  (m) 

Year  of 
optu aiinn 

Mudcl  ice 

CANADA 

NRC/TMD 

Institute  for  Marine  Dynamics 

St.  John’s,  Newfoundland 

94.0  x  12.0x3.0 

1985 

EG/AD 

NRC/Hydraulics  Laboratory 
Ottawa,  Ontario 

18.0  x  7  0x1.2 

1980 

EG/AD/S 

Arctec  Canada  Ltd. 

Calgary,  Alberta 

30.5x7.3x1.4 

1981 

Synthetic 

Arctec  Canada  Ltd. 
Kanata,  Ontario 


30.5  x  4.9  x  1.5 


1977 


Saline 


Table  1  (cont’d).  Existing  ice  testing  facilities.  The  basin  length  given  may 
include  the  length  of  the  trim  tank. 


Main  dimensions 

L  x  B  x  D  (m) 

Year  of 
operation 

Model  ice 

FINLAND 

Wartsila  Arctic 

Research  Center 

Helsinki 

77.3  x  6.5  x  2.3 

1983 

FG-saline 

GERMANY 

HSVA/Hamburgische 

Schiflbau  Versuch. 

Hamburg 

78.0x10.0x2.5/5.0 

30.0x6.0x1.2 

1984 

1971 

Urea-doped 

Saline/Urea 

JAPAN 

SRI/Ship  Research  Institute 

Tokyo 

35.0  x  6.0  x  2.1 

1981 

Saline 

NKK  Tsu  Research  Laboratory 

Tsu  City,  Mie 

20.0x6.0x1.8 

1982 

Urea-doped 

IHI/Ishikawajima- 
Harima  Heavy  Industries 
Yokohama 

7.0  x  3.4  x  0.9 

1982 

Urea-doped 

MSB/Mitsui  Ship 

Building  Company 

Ichihara  City 

5.0  x  1.5  x  0.8 

1985 

Urea-doped 

MH  I/Mitsubishi 

Heavy  Industries 

Nagasaki  City 

20.0  x  9.0  x  2.3 

1986 

Urea-doped 

USA 

ARCTEC  Offshore  Corporation 
Columbia,  Maryland 

25.6x10.7x2.4 

30.5x3.7x1.8 

1982 

1974 

Synthetic 

Saline 

Cold  Regions  Research  and 
Engineering  Laboratory  (CRREL) 
Hanover,  New  Hampshire 

36.5  x  9.1  x  2.4 

1978 

Urea-doped 

Iowa  Institute  of 

Hydraulic  Research 

Iowa  City,  Iowa 

20.0x5.0x1.2 

1980 

Urea-doped 

USSR 

Krylov  Research  Institute 

30.0x6.0x1.8 

1986 

Saline 

Arctic  &  Antarctic 
Research  Institute 
Leningrad 


13.5x1  9x1.3 


1955 


Salir.e 


and  3)  to  investigate  new  concepts  and  de¬ 
signs.  As  a  loose  analogy  with  ice-free  (clear 
water)  tests  either  in  calm  water  or  in  waves, 
model  tests  ’a  ice  can  be  divided  into  two  cate¬ 
gories,  namely  tests  in  level  ice  and  tests  in 
ridges.  At  this  time  procedures  and  tech¬ 
niques  for  tests  in  level  ice  are  fairly  well  es¬ 
tablished  and  documented  in  the  literature 
(e.g.,  Schwarz  1983;  Free,  in  press),  while 
model  tests  in  ridges  are  far  more  recent  and 
their  procedures  are  still  being  developed. 
Model  tests  in  level  ice  have  also  been  the  pri¬ 
mary  concern  of  the  Committee  on  Perfor¬ 
mance  in  Ice-Covered  Waters  of  the  Interna¬ 
tional  Towing  Tank  Conference  (ITTC), 
which  has  only  recently  begun  to  tackle  the 
difficulties  involved  in  model  testing  of  ships 
in  ridges.  The  present  paper  addresses  only 
model  testing  in  level  ice,  where  there  still 
remains  some  controversy,  as  discussed  in 
the  main  body  of  the  paper. 

The  views  expressed  in  this  paper  are  those 
currently  held  by  this  writer.  It  is  likely  not 
only  that  they  will  change  with  time  as  more 
experience  is  gained  and  progress  is  made, 
but  also  that  some  of  them  will  be  challenged 
by  other  researchers  in  this  exciting  field  of 
ice  model  testing. 

GENERAL  CONSIDERATIONS 

As  in  ice-free  model  tests,  the  dependent 
variables  to  be  measured  in  ice  model  tests 
are  the  total  resistance  ( R )  for  EHP  (effective 
horsepower)  tests,  or,  in  propulsion  or  SHP 
(shaft  horsepower)  tests,  the  propeller  thrust 
(Th)  and  propeller  torque  «?),  from  which  oth¬ 
er  ship  propulsion  characteristics  such  as 
shaft  horsepower,  propulsion  efficiency  and 
thrust  deduction  coefficient  can  be  calculat¬ 
ed. 

The  difference  between  model  tests  in  ice- 
free  waters  and  in  level  ice  is  the  presence  of 
a  solid,  flexible,  floating  boundary  at  the  free 
surface.  This  adds  considerably  to  the  diffi¬ 
culty  of  the  tests  by  greatly  increasing  the 
number  of  independent  variables  or  parame¬ 
ters  that  need  to  be  scaled.  These  variables, 
listed  in  Table  2,  include  not  only  the  usual 
ones,  namely  ship  geometric  characteristics, 
ship  speed,  propeller  speed,  water  properties 
(density  and  viscosity)  and  acceleration  of 
gravity,  but  also  the  ice  characteristics, 


namely  thickness,  density  and  mechanical 
properties,  and  last  but  not  least  the  ice-hull 
friction  coefficient.  All  of  these  parameters 
must  be  scaled  to  satisfy  similitude.  There¬ 
fore,  ship  modeling  in  ice  involves  hydrody¬ 
namics  as  well  as  mechanics  of  materials  to 
account  for  ice  deformation  and  rupture.  It  is 
thus  strongly  dependent  on  the  progress  made 
in  the  field  of  ice  mechanics.  However,  the 
term  “ice”  encompasses  many  different  ma¬ 
terials:  freshwater  ice,  first-year  sea  ice  and 
multi-year  sea  ice,  to  name  a  few,  with  vary¬ 
ing  crystalline  structure  and  corresponding 
behavior  under  load.  In  spite  of  the  difficul¬ 
ties  outlined  above,  the  art  of  modeling  ship- 
ice  interaction  has  made  significant  advanc¬ 
es  in  the  past  decade  to  the  point  that  results 
from  model  tests  can  greatly  aid  the  naval 
architect  in  the  design  of  an  ice-transiting 
ship. 

Perfect  physical  modeling  of  any  phenom¬ 
enon  requires  exact  geometric,  kinematic 
and  dynamic  similitude.  Exact  similitude  is 
ensured  when  all  independent,  dimension¬ 
less  parameters  involved  take  the  same  val¬ 
ues  at  model  scale  as  at  full  scale.  For  the 
particular  problem  of  sbip-ice  interaction, 
these  dimensionless  parameters  are  listed  in 
Table  2,  where  the  selected  nondimensional- 
izing  variables  are  ice  thickness  (h{),  water 
density  (p)  and  gravity  (g). 

Since  the  gravity,  g,  cannot  usually  be  al¬ 
tered  in  the  model  tests,  and  the  water  in  the 
tank  has  nearly  the  same  properties  (density 
and  viscosity)  as  sea  water,  it  is  well  known 
that  identity  of  the  Froude  number  and  the 
Reynolds  number  between  model  and  proto¬ 
type  cannot  be  achieved.  Since  tests  in  ice  are 
usually  made  at  relatively  low  velocities, 
corresponding  to  a  maximum  of  6  to  8  knots  (3 
to  4  m/s)  at  full  scale,  the  resistance  compo¬ 
nent  due  to  the  water  itself  is  usually  a  small 
percentage  of  the  total  resistance  and  either 
can  be  neglected  altogether,  or  accommodated 
with  a  minor  correction  of  the  test  results.  In 
this  respect,  it  should  be  noted  that  the  pres¬ 
ence  of  turbulence  stimulators  at  the  bow,  a 
common  practice  in  open  water  tank  tests,  is 
not  acceptable  in  ice  tests  since  even  a  small 
protrusion  in  the  model  hull  will  unrealisti¬ 
cally  increase  the  ice  resistance.  The  main 
concern  of  the  engineer  and  scientist  in 
charge  of  ice  tests  is,  therefore,  to  model  as 
correctly  as  possible  the  ice  properties  accord¬ 
ing  to  the  modeling  law  selected. 
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Table  2.  Variables  governing  ship-ice  interaction. 


Dimensional 

Dimensionless 

Type  of  variables 

variables 

parameter 

Ship 

Length 

L 

LI  h  i 

characteristics 

Beam 

B 

B/h, 

Draft 

D 

D/hi 

Propeller  diameter 

d 

dJk  i 

Speed 

V 

Fn  =  V  /  ‘fgh]  or 

Fd  =  V  /V  (Ap  I p)gkt 

Propeller  speed 

n 

J=VI  nd 

Gravity 

g 

1 

Water 

Density 

p 

1 

properties 

Viscosity 

p 

pVL 

P 

Thickness 

1 

Density 

Pi 

Pit? 

Poisson’s  ratio 

V 

u 

Modulus 

E 

E/pghi  or 

El  Of  or 
c  n  =  PV2/E 

Ice 

Characteristic  length 

K 

Ic/hi 

properties 

Flexural  strength 

°r 

O  [Ipghi 

Comp,  strength 

CTC 

CTc/or 

Shear  strength 

Ob 

OJOf 

Fracture  toughness 

*lc 

In  =PiV  <IJKU 

15 

or  KlcPighi 

Ice-hull  friction  factor 

fi 

fi 

NOTES: 

The  ice  characteristic  length,  lc,  is  defined  by: 


1  = 


Eh . 


12  pa  (i  -  v  ) 


1/4 


Fn  is  the  Froude  number  based  on  ice  thickness. 
Fa  is  the  densimetric  Froude  number. 

J  is  the  advance  coefficient. 

Cn  is  the  Cauchy  number. 

/„  is  the  “Ice”  number  propose,  by  Atkins  (1975). 


The  current  practice  in  ice  modeling  is  to 
follow  the  Froude  scaling  law.  If  the  geomet¬ 
ric  scale  is  selected  first,  X  =  Lp/Lm  =  proto¬ 
type  length/model  length,  the  ratio  of  proto¬ 
type  velocity  to  model  velocity  is  given  by 

*v  =  VVm=  X°5- 

Conversely,  it  can  be  envisioned  that  the 
available  towing  capabilities  will  limit  the 
range  of  model  speeds  and  therefore  dictate 
the  velocity  scale  Xv.  The  geometric  scale  of 
the  model  would  then  be  selected  as  X  =  X^. 
Other  factors  that  also  affect  the  choice  of  the 
model  geometric  scale  are:  model  tank  size, 
model  cost,  minimum  model  ice  strength  for 
which  test  results  are  reliable,  characteristics 
of  available  dynamometers  for  measuring 
resistance,  thrust  and  torque,  or  even  size  of 
available  model  propellers.  The  model  scale 
finally  selected  is  a  compromise  between  the 
often  conflicting  above  considerations. 

Once  the  geometric  scale  has  been  selected, 
the  scales  of  all  the  other  test  variables  and 
thus  their  range  of  values  to  be  achieved  in  the 
model  tests  are  specified  by  the  requirements 
that  the  dimensionless  parameters  in  Table  2 
take  the  same  values  in  the  model  as  in  the 
prototype.  This  implies  in  particular  that  ice 
thickness,  ice  modulus,  and  ice  strength  be 
scaled  by  the  geometric  scale  X,  while  its  frac¬ 
ture  toughness  should  be  scaled  by  X1-5.  The 
parameters  v  (ice  Poisson’s  ratio)  and  f j 
(hull-ice  friction  coefficient),  which  have  no 
dimensions,  should  have  the  same  values  in 
the  model  as  in  the  prototype.  The  first  re¬ 
quirement  is  to  know  the  range  of  prototype 
values  for  the  ice  properties  that  must  be  mod¬ 
eled.  A  brief  discussion  on  ice  properties  is 
given  in  the  following  section. 

ICE  PROPERTIES 

Much  research  has  been  devoted  to  the  field 
and  laboratory  measurements  of  ice  mechan¬ 
ical  properties.  Relevant  information  can  be 
found  in  Weeks  and  Assur  (1967),  Lavrov 
(1969),  Schwarz  and  Weeks  (1977),  Tryde 
(1979),  Mellor  (1983),  in  the  chapters  on  ice 
physics  and  ice  mechanics  of  books  devoted  to 
ice  engineering  (Michel  1978,  Ashton  1986), 
and  in  specialized  scientific  journals  such  as 
Journal  of  Glaciology,  Cold  Regions  Science 


and  Technology ,  and  others.  The  proceed¬ 
ings  of  specialty  conferences  such  as  the  Ice 
Symposia  of  the  International  Association  for 
Hydraulic  Research  (IAHR),  Port  and  Ocean 
Engineering  under  Arctic  Conditions 
(POAC)  conferences,  and  Offshore  Mechan¬ 
ics  and  Arctic  Engineering  (OMAE)  confer¬ 
ences  are  an  important  source  of  up  to  date  in¬ 
formation. 

Even  a  cursory  survey  of  the  literature 
shows  that  ice  mechanical  properties  are 
highly  variable.  The  values  measured  de¬ 
pend  upon  the  type  of  ice  (i.e.,  freshwater  ice 
or  saline  ice)  and  its  structural  characteris¬ 
tics  (i.e.,  grain  size,  crystallography  and  po¬ 
rosity  or  brine  volume  in  the  case  of  sea  ice), 
temperature,  rate  of  strain,  and  even  on  the 
measurement  techniques.  To  address  this 
latter  problem,  the  ITTC  Committee  on  Per¬ 
formance  in  Ice-covered  Waters  and  the 
IAHR  Committee  on  Ice  Problems  are  at¬ 
tempting  to  standardize  the  experimental 
techniques  for  measuring  ice  properties  both 
in  the  field  and  in  the  laboratory  (Schwarz 
1979).  It  should  be  mentioned  at  the  outset  that 
measurement  techniques  do  vary  between 
field  and  laboratory  conditions  (Timco 
1981b),  and  that  the  values  attributed  to  the 
natural  and  even  more  so  to  the  model  ice 
should  be  treated  more  as  indices  than  as  ac¬ 
tual  mechanical  properties  of  ice  as  a  materi¬ 
al. 

Modulus,  E 

The  ice  modulus,  E,  used  in  engineering  is 
actually  the  initial  tangent  or  effective  modu¬ 
lus  determined  from  the  load-displacement 
curve  of  a  plate,  a  cantilever  beam  or  a  sim¬ 
ply  supported  beam.  The  rate  of  load  applica¬ 
tion  should  be  high  enough  for  creep  deforma¬ 
tion  of  the  ice  to  be  negligible  or  corrections 
need  to  be  introduced  when  interpreting  the 
test  data.  This  effective  modulus  is  known  to 
vary  with  ice  salinity  and  temperature,  i.e., 
with  porosity.  Vaudrey  (1977)  gives  the  fol¬ 
lowing  relationship  between  E  and  porosity, 
p,  of  sea  or  saline  ice: 

E  (GPa)  =  5.25  - 13  Vp  (p  <  0.1). 

For  freshwater  ice  (p= 0),  the  effective  mod¬ 
ulus  predicted  by  the  above  equation  is  5.25 
GPa,  which  is  within  the  measured  range  of  4 
to  6  GPa.  The  modulus  can  be  as  low  as  1  GPa 


for  warm  sea  ice  with  a  porosity  of  0.1  or 
greater  (brine  volume  of  100%oo  or  more). 
Other  reseachers  claim  a  linear  variation  of 
E  withp. 

The  ice  Young’s  or  elastic  modulus  is 
measured  usually  by  high  frequency  vibra¬ 
tional  methods,  i.e.,  by  measuring  the  rate  of 
propagation  of  small-amplitude,  high- 
frequency  waves  or  pulses  in  the  ice.  The  val¬ 
ues  obtained  in  this  manner  are  nearly  twice 
those  of  the  effective  modulus. 


Characteristic  length,  lc 

When  analyzing  plate  or  beam  test  results 
for  the  determination  of  E,  an  often  implicit 
calculation  is  that  of  the  characteristic 
length,  lc.  This  length  is  a  measure  of  the  size 
of  the  deformation  zone  of  the  ice  plate  or 
beam  when  subjected  to  a  vertical  static  load. 
The  two  quantities  E  and  lc  are  related  by 


for  a  plate,  and 


for  a  beam. 

To  determine  E  from  lc,  one  has  to  take  the 
ice  thickness,  hit  into  account.  Also,  it  can  be 
seen  that  an  error  of  18%  in  Zc  leads  to  a  100% 
error  in  E.  Therefore,  scaling  the  character¬ 
istic  length  may  be  preferable  to  scaling  E.  In 
measuring  lc  by  the  plate  method  (Sodhi  et  al. 
1982)  for  both  lake  ice  and  first-year  sea  ice, 
Sodhi  et  al.  (1985)  found  that  the  ratio  IJh^  for 
lake  ice  was  on  the  order  20,  while  for  sea  ice 
it  was  in  the  range  of  10  to  15. 


Flexural  strength,  <7t 

From  results  of  cantilever  beam  tests  and 
simply  supported  beam  tests,  it  was  found  that 
Of  decreased  from  approximately  1  GPa  for 
freshwater  ice  to  about  300  kPa  for  sea  ice  with 
a  porosity  of  0.2  or  greater.  Thus  the  ratio  of 
effective  modulus  to  flexural  strength,  E/a(, 
varies  approximately  between  5000  (for  fresh¬ 
water  ice  and  cold  first-year  sea  ice)  to  2000 
(for  warm,  weak  sea  ice).  Vaudrey  (1977),  in 
agreement  with  Weeks  and  Assur  (1967),  pro¬ 
posed  the  following  relationship  for  flexural 


strength  of  columnar  S2  sea  ice  in  terms  of 
the  ice  porosity,  p: 

CTf=  960(1  -  2  Vp)kPa 

valid  forp  <0.11.  Forp  >0.11,  crf  remains  ap¬ 
proximately  constant,  equal  to  250-350  kPa. 

For  columnar  freshwater  ice,  it  was  found 
that  the  flexural  strength  obtained  from  beam 
tests  (cantilever  or  simply  supported  beams) 
depended  on  the  direction  (upward  or  down¬ 
ward)  of  load  application  (Gow  et  al.  1978). 
The  values  of  Of  obtained  with  the  ice  top  fiber 
in  tension  is  1.25  to  2  times  that  measured 
with  the  bottom  fiber  in  tension.  This  differ¬ 
ence  is  attributed  to  the  anisotropic  crystal¬ 
line  structure  of  columnar  ice,  where  the 
grain  size  increases  from  top  to  bottom.  For 
sea  ice,  however,  there  is  no  apparent  differ¬ 
ence  in  strength  between  the  two  loading  di¬ 
rections  (Tabata  et  al.  1967,  Kayo  et  al.  1983). 
Also,  the  flexural  strength  of  freshwater  ice 
depends  on  the  measurement  technique,  sim¬ 
ply  supported  beam  tests  yielding  values  as 
much  as  twice  those  obtained  from  cantilever 
beam  tests  (Gow  et  al.  1978,  Timco  and  Fre- 
derking  1982,  Timco  1985).  No  such  effect  of 
measurement  technique  on  flexural  strength 
was  observed  for  sea  ice. 

In  view  of  the  above,  it  is  worth  emphasiz¬ 
ing  again  that  should  be  considered  as  a 
flexural  strength  index  and  that,  when 
presenting  data  on  crr,  one  should  also  indi¬ 
cate  the  measurement  methods  and  condi¬ 
tions. 

Compressive  strength,  ac 

Since  the  maximum  ice  forces  exerted  on  a 
structure  occur  when  ice  fails  in  crushing, 
the  compressive  strength  is  the  property  that 
has  received  the  most  attention  from  re¬ 
searchers.  Through  careful  testing,  it  has 
been  determined  that  this  property  depends  on 
the  following  parameters:  ice  type,  grain 
size,  temperature  or  porosity  (or  both),  strain 
rate  and  test  conditions  (i.e.,  uniaxial  test  vs 
triaxial  test,  relative  orientation  of  force  ap¬ 
plication  to  crystalline  structure,  etc.)  (Tim¬ 
co  and  Frederking  1986).  For  the  purpose  of 
the  present  paper,  i.e.,  for  modeling  of  ships 
in  ice  where  the  loading  rate  is  high,  it  may 
be  sufficient  to  say  that  crc  in  the  horizontal  di¬ 
rection  is  3  to  5  times  <rf,  while  in  the  vertical 
direction  it  can  be  as  much  as  10  times  <rf. 
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Shear  strength,  crs 

Data  on  ice  shear  strength  are  very  few  pri¬ 
marily  because  no  fully  reliable  methodology 
to  measure  this  ice  property  has  yet  been  de¬ 
veloped.  On  the  basis  of  limited  field  experi¬ 
ments,  Frederking  and  Timco  (1986)  have 
proposed  the  following  equation  relating  the 
vertical  shear  strength  of  columnar  ice  to  ice 
porosity,  p: 

<rs(kPa)  =  1500  [1  - 1.60^] 

which  is  in  agreement  with  the  commonly  ac¬ 
cepted  assumption  that  <rs  is  1.5  to  2  times  the 
ice  flexural  strength. 

Poisson’s  ratio,  v 

To  quote  Mellor  (1983,  p.  53)  “the  value  of  v 
for  non-saline  ice  of  very  low  porosity  is 
about  0.33  ±  0.03,  and  variation  with  porosity 
in  sea  ice  is  likely  to  be  within  the  limits  of 
uncertainty  for  the  pure  ice  value.”  It  is  com¬ 
mon  to  take  a  value  of  u  =  0.3  for  sea  ice  and 
model  doped  ice  (either  saline  or  urea  ice). 

Fracture  toughness,  KIc 

Until  very  recently,  fracture  toughness 
had  been  a  property  totally  ignored  in  the 
analysis  of  ship-ice  interaction.  It  is,  howev¬ 
er,  receiving  increasing  attention  and  its 
importance  in  ship  modeling  in  ice  is  being 
investigated  (Jones  1986,  Parsons  et  al.  1986). 
It  characterizes  a  material’s  resistance  to 
crack  growth.  It  is  likely,  therefore,  to  control 
the  size  of  the  resulting  pieces  when  the  mate¬ 
rial  fails  in  fracture  (fragmentation).  De¬ 
fects  or  flaws  in  the  material  will  decrease  its 
fracture  toughness.  Air  bubbles  in  freshwater 
ice  and  brine  pockets  in  sea  or  doped  ice  are 
such  defects.  This  material  property  is  possi¬ 
bly  the  most  difficult  to  measure  experimen¬ 
tally.  It  is  therefore  not  surprising  that  the 
data  published  in  the  literature  for  freshwater 
and  sea  ice  exhibit  considerable  scatter  (see 
Mellor  [1983]  for  references).  The  data  avail¬ 
able  in  the  literature  show  values  for  the  frac¬ 
ture  toughness  that  vary  for  freshwater  ice 
from  60  to  300  kN  m'3/2,  and  for  sea  ice  from 
30  to  140  kN  nr 3/2 

Density,  pi 

While  the  density  of  pure  ice  has  long  been 
established  at  917  kg/m3,  that  of  field  and  la¬ 


boratory  ice  is  affected  by  the  presence  of  im¬ 
purities  such  as  air  bubbles  or  brine  pockets. 
The  density  of  bubbly  freshwater  ice  can  be  as 
low  as  890  kg/m3.  For  sea  ice,  the  usual  range 
of  density  is  from  about  890  kg/m3  for  cold, 
low  salinity  ice  to  930  kg'm3  for  ice  with  high 
brine  content. 

Ice-hull  friction  factor, 

When  a  ship  penetrates  into  ice,  a  signifi¬ 
cant  component  of  the  resistance  to  motion 
may  come  from  e  friction  of  the  ice  on  the 
hull.  This  frictn  ial  force  is  assumed  to  be 
proportional  to  the  force  exerted  by  the  ice  nor¬ 
mal  to  the  ship-hull,  and  the  coefficient  of 
proportionality  is  the  friction  coefficient  /";. 
Numerous  studies  (e.g.,  Ryvlin  1973,  Cala¬ 
brese  et  al.  1980,  Oksanen  1983,  Saeki  et  al. 
1984,  Akkok  et  al.  1986,  Tatinclaux  et  al. 
1986)  have  been  conducted  to  determine  the  ef¬ 
fect  of  various  parameters  on  fx.  It  was  found 
that  ice  hardness,  relative  velocity  between 
ice  and  hull,  wetting  properties  of  the  hull, 
etc.,  could  affect  fx,  in  addition  to  hull  rough¬ 
ness  characteristics  such  as  roughness  aver¬ 
age  and  morphology  of  the  roughness  ele¬ 
ments.  Here  again,  the  methodology  and  ap¬ 
paratus  used  to  measure  fx  may  affect  the  val¬ 
ues  obtained.  Most  studies  were  made  in  the 
laboratory  or  under  controlled  conditions. 
The  few  direct  measurements  of  fx  on  real 
ships  (Hoffman  1985,  Liukkonen  1985,  Enk- 
vist  and  Mustamaki  1986)  gave  inconclusive 
and  even  inconsistent  results.  Low  friction 
coatings,  the  best  known  of  which  is  Inerta- 
160,  have  been  and  are  being  developed.  In 
practice,  it  is  generally  assumed  that  a  new 
ship  hull  has  a  friction  coefficient  of  about 
0.1. 

Incomplete  as  it  may  be,  the  above  discus¬ 
sion  should  give  a  fair  idea  of  the  difficulty 
in  determining  the  ice  properties.  Values  for 
these  properties  should  always  be  accompa¬ 
nied  by  a  careful  description  of  the  test  meth¬ 
ods  and  conditions  used  to  obtain  them,  as 
well  as  information  on  ice  structure,  salini¬ 
ty,  etc.  It  is  worth  repeating  that,  in  fact,  the 
values  measured  do  not  actually  represent  the 
material  properties  of  the  particular  ice  test¬ 
ed,  but  should  be  considered  as  indices  to  be 
used  for  modeling  or  for  comparison  between 
model  test  results. 
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MODEL  ICE 


As  discussed  previously,  once  the  geomet¬ 
ric  scale  of  the  model  has  been  selected,  the 
mechanical  properties  of  the  model  ice  ought 
to  be  in  a  specified  ratio  with  those  of  the  field 
ice  to  be  modeled.  There  are  currently  five 
types  of  model  ice  used  in  the  various  tanks, 
as  indicated  in  Table  1,  namely  synthetic  ice, 
saline  columnar  ice,  columnar  carbamide- 
doped  (urea)  ice,  fine-grained  ice  and  EG/ 
AD/S  ice.  Except  for  the  first,  all  types  of  mod¬ 
el  ice  are  grown  from  water  baths  containing 
one  or  more  chemicals  or  dopants.  These  do¬ 
pants  are  partially  trapped  between  the  ice 
grains,  thereby  reducing  the  strength  of  the 
ice.  Further  adjustment  in  the  strength  char¬ 
acteristics  of  the  doped  ice  can  be  made  some¬ 
what  by  varying  the  growth  temperature  and 
primarily  by  raising  the  ice  sheet  tempera¬ 
ture  to  within  a  few  degrees  of  the  melting 
temperature  for  a  certain  amount  of  time 
(warmup  period)  prior  to  the  tests.  An  excel¬ 
lent  discussion  of  the  morphological  proper¬ 
ties  of  ice  grown  from  a  doped  solution  was 
given  by  Timco  (1979).  This  section  discuss¬ 
es  the  main  properties  of  these  types  of  model 
ice  with  their  advantages  and  disadvantages 
as  can  be  gathered  from  the  open  literature. 

Synthetic  ice 

The  more  widely  known  synthetic  ice 
(MOD  ice)  is  a  wax-based  mixture  prepared 
in  a  liquid  form  and  poured  on  the  surface  of 
a  traditional  towing  tank  where  it  is  allowed 
to  set.  The  exact  composition  of  the  mixture  is 
adjusted  to  achieve  as  closely  as  possible  the 
required  density  and  mechanical  properties. 
Such  a  synthetic  ice  was  used  as  early  as  1970 
(Crago  et  al.  1970).  Its  rights  were  acquired  by 
Arctec  Inc.  (now  Arctec  Offshore  Corpora¬ 
tion)  of  Columbia,  Maryland,  where  it  has 
been  further  refined.  It  is  a  combination  of  a 
“secret  ingredient”  and  flat,  round  plastic 
beads  or  miniature  disks,  3  to  4  mm  across. 
Since  the  exact  composition  of  this  synthetic 
ice  is  proprietary,  no  other  ice  tank  besides 
those  of  Arctec  Offshore  Corporation  in  the 
U.S.A.  and  of  Arctec  Canada,  Ltd.,  has  any 
experience  with  this  model  ice  and  there  has 
been  no  independent  study  of  its  properties. 
Schultz  and  Free  (1984)  presented  the  results 
of  measuring  the  properties  of  40  model  ice 
sheets.  The  ice  thicknesses  ranged  from  10  to 


60  mm,  with  most  ice  sheets  (33  out  of  40  or 
83%)  between  15  and  35  mm  thick.  The 
flexural  strength  ranged  from  as  low  as  10 
kPa  to  about  80  kPa;  however,  34  out  of  40  ice 
sheets  (85%)  had  a  flexural  strength  between 
10  and  20  kPa.  The  range  of  £/<rf  cited  is  700  to 
5000  with  an  average  of  1450  and  with  90%  of 
the  data  points  lying  between  700  and  2000. 
The  ratio  of  compressive  strength  to  flexural 
strength,  ac/af,  ranged  from  0.6  to  2.3  with 
an  average  of  1.14,  which  is  low  compared  to 
that  of  sea  or  freshwater  ice.  Finally,  the  ice- 
model  hull  friction  factor,  was  given  as 
ranging  from  0.07  to  0.28,  which  covers  the 
entire  realistic  full-scale  range  that  can  be 
expected.  Schultz  and  Free  did  not  give  val¬ 
ues  for  the  density,  the  shear  strength  or  the 
fracture  toughness  of  the  synthetic  model  ice. 
They  did,  however,  mention  that  “the  size  of 
the  broken  ice  pieces  resulting  from  ship  or 
structure  interaction  with  the  synthetic  ice  ... 
appears  to  compare  very  favorably  with  the 
corresponding  piece  size  observed  in  field 
trials  of  icebreaking  vessels  and  structures 
(p.  239).”  They  also  cite  one  definite  advan¬ 
tage  of  the  synthetic  ice  over  doped  model  ice: 
it  is  handled  at  room  temperature  and  its 
properties  remain  unchanged  for  days. 

Columnar  saline  ice 

The  first  dopant  used  to  create  model  ice  in 
a  refrigerated  ice  tank  was  sodium  chloride. 
By  varying  the  salt  concentration  in  the  wa¬ 
ter  bath  and  consequently  the  brine  content 
within  the  ice,  almost  any  required  reduction 
in  flexural  strength  can  be  achieved  at 
growth  temperature,  which  is  normally  well 
below  0°C.  This  is  still  the  technique  used  in 
the  ice  tanks  of  the  Arctic  and  Antarctic  Re¬ 
search  Institute  in  Leningrad,  U.S.S.R. 
However,  high  salinity  model  ice  has  an  ex¬ 
tremely  low  strain  modulus  E,  and  the  ratio 
E/Of  is  usually  no  greater  than  500.  This  ice 
is  therefore  excessively  plastic,  and  breaks 
into  floes  that  are  much  larger  relative  to 
those  observed  in  the  field.  Schwarz  (1975, 
1977)  suggested  use  of  a  low  salinity  ice 
grown  from  a  water  bath  with  a  salt  concen¬ 
tration  of  1.6  %c  or  less.  The  method  suggest¬ 
ed  involves  three  steps: 

1.  When  the  water  bath  has  nearly  reached 
its  freezing  temperature,  a  fine  mist  is 
sprayed  into  the  cold  air  above  the  tank.  The 
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mist  freezes  in  fine  crystals  that  settle  on  the 
water  surface  to  initiate  the  ice  sheet.  This 
wet'Seeding  technique  results  in  a  columnar 
ice  sheet  with  small  crystals,  one  that  closely 
resembles  the  crystalline  structure  of  sea  ice. 

2.  The  ice  is  grown  at  as  low  a  temperature 
as  possible  (on  the  order  of  -20°C),  primarily 
to  reduce  the  growth  time  to  a  minimum,  but 
also  to  minimize  the  thickness  of  the  nearly 
freshwater  top  layer. 

3.  Once  the  desired  ice  thickness  is  nearly 
reached,  the  temperature  in  the  ice  tank  is 
raised  to  about  0°C.  During  this  warmup  peri¬ 
od,  the  brine  pockets  in  the  ice  melt,  thereby 
increasing  the  liquid  volume  within  the  ice 
and  decreasing  the  flexural  strength.  With 
this  procedure,  Schwarz  showed  that  values  of 
the  ratio  E/a f  close  to  2000  were  possible  as 
long  as  the  ice  strength  was  above  40  to  50 
kPa.  These  were  confirmed  in  an  indepen¬ 
dent  study  by  Timco  (1981a).  However,  when 
further  reduction  in  strength  was  required, 
i.e.,  for  a  geometric  scale  greater  than  15  to 
10,  the  ratio  E/af  could  become  as  low  as  500. 
Another  disadvantage  of  salt  as  a  dopant  is 
that  it  is  highly  corrosive.  Its  use  leads  to  high 
maintenance  costs  for  the  steel  components  of 
the  ice  tank,  towing  carriage,  and  other 
equipment. 

Columnar  carbamide  ice 

After  an  extensive  study  of  possible  do¬ 
pants,  Timco  (1979,  1980,  1981a)  suggested 
carbamide  (or  urea)  as  an  alternative  to  sodi¬ 
um  chloride  for  high  modeling  scale  (>20). 
From  test  results  made  on  4-  to  5-cm  thick  ice 
sheets  grown  from  a  1.3%  by  weight  urea  so¬ 
lution  in  water,  Timco  concluded  that  carba¬ 
mide  was  superior  to  low-salinity  saline  ice 
and  that  scale  factors  up  to  40  were  possible 
with  this  new  doped  ice.  Over  the  range  of 
flexural  strength  from  15  to  70  kPa,  the  ratio 
E/a{  remained  larger  than  2000,  the  compres¬ 
sive  strength  in  the  vertical  direction  varied 
between  110  and  250  kPa  and  that  in  the  hori¬ 
zontal  direction  varied  between  110  and  160 
kPa.  The  shear  strength  under  vertical  load¬ 
ing  increased  from  30  kPa  at  a{=  15  kPa  to  70 
kPa  at  <Tf  =  70  kPa,  while  that  under  horizon¬ 
tal  loading  increased  from  35  to  65  kPa.  The 
density  of  carbamide  ice  was  found  to  be  on 
the  order  of  0.95  ©'em3.  Measurements  made 


at  CRREL  yielded  values  of  0.93  g/cm3.  A 
significant  advantage  of  carbamide  over  salt 
is  that  it  is  not  corrosive.  It  is  also  a  very  com¬ 
mon,  fairly  inexpensive  chemical  used,  in 
particular,  as  a  fertilizer. 

Many  ice  tanks  —  Arctec  Offshore  Corpo¬ 
ration,  Arctec  Canada,  CRREL,  HSVA,  NKK 
—  have  adopted  carbamide  ice  as  model  ice. 
Hirayama  (1983),  in  an  extensive  study  of 
carbamide  ice  sheets  grown  from  a  0.95% 
urea  solution  in  water,  confirmed  most  of 
Timco’s  findings.  However,  in  the  experi¬ 
ence  of  this  writer,  while  values  of  E/a f  of  2000 
are  possible  for  relatively  thick  and  strong 
ice,  values  of  1000  to  1500  are  more  common. 
For  thin,  weak  ice,  i.e.,  for  ice  sheets  less 
than  3  cm  in  thickness  and  with  a  flexural 
strength  less  than  30  kPa,  the  ratio  E/<t(  may 
even  drop  below  1000.  Correspondingly,  while 
the  ratio  ljht  is  on  the  order  of  12  for  ice  sheets 
4  cm  thick  or  more  and  with  a  flexural 
strength  of  40  kPa  or  above,  it  may  become 
less  than  10  for  thin,  weak  ice  (i.e.,  2  cm  thick 
and  with  ar  =  20  kPa).  For  these  reasons,  this 
writer  discourages  users  of  the  CRREL  ice 
tank  from  conducting  tests  at  thicknesses  of 
less  than  2  cm  or  ice  flexural  strengths  less 
than  20  kPa.  In  this  writer’s  experience,  the 
results  of  such  tests  would  be  unreliable. 

Another  disadvantage  of  urea-doped  ice  is 
its  two-layer  structure:  the  top,  incubation 
layer  can  be  3  to  10  mm  in  thickness  and  is 
significantly  stronger  than  the  bottom,  den¬ 
dritic  layer.  For  that  reason,  the  flexural 
strength  measured  with  the  load  applied  up¬ 
ward  (bottom  surface  in  tension)  is  only  about 
40%  of  that  measured  with  the  load  applied 
downward  (top  surface  in  tension).  The 
thickness  of  the  top  layer  can  be  minimized 
by  growing  the  ice  at  the  lowest  room  tempera¬ 
ture  possible  and  ensuring  that  the  water  tem¬ 
perature  at  seeding  is  nearly  equal  to  its 
freeezing  point. 

Finally,  Timco  (1985)  measured  values  of 
the  critical  stress  intensity  factor  7fIc  from  an 
average  of  about  6  kPa  m1/2  at  <xr  =  40  kPa,  to 
an  average  of  17  kPa  m1/2  at  <rf  =  130  kPa. 
When  extrapolated  to  full  scale  according  to 
flexural  strength  ratio,  these  values  of  Jfjc  are 
significantly  greater  than  the  range  of  30  to 
100  kPa  m1/2  for  sea  ice  with  a  porosity  of  0.13 
and  that  of  90  to  150  kPa  m1/2  for  ice  with  a  po¬ 
rosity  of  about  0.02  (Urabe  and  Yoshitake 
1984,  Timco  and  Frederking  1983).  It  can  be 


conjectured  that  the  hard  top  layer  of  carba¬ 
mide  ice  may  be  at  least  partially  responsible 
for  its  high  critical  stress  intensity  factor. 
Another  reason  might  be  that  model  ice  has 
fewer  flaws  than  real  ice.  This  high  fracture 
toughness  of  carbamide  ice,  even  with  a  high 
ratio  of  £/<rf,  may  be  one  reason  that  the  size  of 
the  ice  floes  produced  with  this  model  ice  re¬ 
main  relatively  larger  than  those  observed 
during  full-scale  trials. 

While  in  the  opinion  of  this  writer  carba¬ 
mide  ice  was  a  definite  improvement  over  sa¬ 
line  ice,  it  is  not  yet  the  ideal  material  for 
modeling  ice-structure  interactions. 

Fine-grained  ice 

Almost  concurrently  with  the  development 
of  carbamide  ice  at  the  National  Research 
Council  in  Ottawa,  Enkvist  and  co-workers 
from  Wartsila  Arctic  Research  Center 
(WARC)  in  Finland  were  developing  anoth¬ 
er  model  ice  (Enkvist  and  Makinen  1984).  It 
is  a  saline,  fine-grained  ice  similar  in  struc¬ 
ture  to  snow  ice  (ice  type  T-l)  grown  from  a 
2%  NaCl  solution.  Its  structure  is  nearly  ho¬ 
mogeneous,  as  opposed  to  the  columnar  struc¬ 
ture  of  sea  ice  and  carbamide  ice.  It  is  grown 
by  continuously  spraying  a  mist  of  tank  wa¬ 
ter  above  the  water  surface.  Since  the  air  tem¬ 
perature  is  well  below  0°C,  the  mist  freezes  in 
small  crystals  that  deposit  on  the  water  sur¬ 
face.  This  procedure  requires  that  the  spray 
nozzles  be  mounted  on  a  carriage  that  contin¬ 
uously  travels  back  and  forth  along  the  ice 
tank.  Careful  control  of  the  carriage  speed, 
nozzle  arrangement  and  discharge  is  neces¬ 
sary  to  create  an  ice  sheet  of  uniform  thick¬ 
ness.  Newly  formed  ice  sheets  have  usually 
too  low  a  strength,  and  additional  freezing 
without  spraying  may  be  necessary.  Accord¬ 
ing  to  Enkvist  and  Makinen,  the  ratio  E/a f  ( E 
measured  by  the  plate  deflection  method)  is 
on  the  order  of  1200  ±  200  for  ice  thicknesses  of 
20  to  40  mm  and  flexural  strengths  from  13  to 
28  kPa.  According  to  its  developers,  the  uni¬ 
axial  compressive  strength  of  this  model  ice 
is  somewhat  low. 

According  to  Enkvist  and  Makinen,  the 
main  advantage  of  this  ice,  dubbed  FG-ice,  is 
that  it  is  more  brittle  than  either  saline  or  car¬ 
bamide  ice,  and  therefore  breaks  into  small¬ 
er  floes  than  either.  The  resulting  icebreak¬ 
ing  pattern  and  flow  of  ice  pieces  along  a 
ship’s  hull  are  similar  to  full-scale  observa¬ 


tions.  Also,  the  rate  of  growth  of  this  ice  is  5  to 
7  mm/hour  compared  to  2  to  3  mm/hour  for  the 
other  types  of  model  ice,  which  permits  a 
much  faster  turn-around  time  in  testing. 

Finally,  the  density  and  friction  factor  of 
the  FG-ice  are  said  to  be  comparable  to  that  of 
columnar  sea  ice  or  model  ice.  No  measure¬ 
ments  of  its  shear  strength  or  critical  stress 
intensity  factor  were  reported. 

No  independent  testing  of  this  fine¬ 
grained  model  ice  confirming  the  above  re¬ 
sults  is  known  to  this  writer.  A  similar  fine¬ 
grained  ice  is  now  in  use  at  the  Nippon  Ko- 
kan  K.K.  ice  tank,  where  carbamide  rather 
than  salt  is  the  dopant.  The  growth  procedure 
of  FG-ice  at  NKK  is  also  somewhat  different 
from  that  developed  by  Wartsila’s  research¬ 
ers,  but  results  in  a  comparable  model  ice. 

EG/AD/S  model  ice 

Timco  (1986)  introduced  a  new  type  of  mod¬ 
el  ice  for  use  in  refrigerated  ice  tanks.  The 
ice  is  grown  from  an  aqueous  solution  of 
three  dopants,  ethylene  glycol  (EG),  aliphatic 
(very  low  sudsing)  detergent  (AD)  and  sugar 
(S).  According  to  Timco,  the  ethylene  glycol 
serves  as  a  low  molecular  weight  dopant  that 
is  trapped  in  the  ice  as  impurity  pockets  to  re¬ 
duce  the  strength  of  the  ice.  The  aliphatic  de¬ 
tergent  reduces  the  surface  tension  of  the  so¬ 
lution  and  allows  more  of  the  ethylene  glycol 
to  be  trapped  by  the  ice.  The  resulting  ice 
should  be  nearly  single-layered,  as  opposed  to 
saline  or  carbamide  ice.  Finally,  the  sugar 
inhibits  the  lateral  growth  of  the  ice  platelets 
to  produce  a  fine-grained  ice. 

The  methodology  to  grow  and  temper  this 
new  ice  to  achieve  prescribed  thickness  and 
mechanical  properties  is  basically  the  same 
as  for  carbamide  ice.  Tests  were  carried  out 
by  Timco  (1986)  on  4-cm-thick  ice  sheets 
grown  from  a  solution  with  concentrations  of 
the  three  chemicals  of  0.46/0.032/0.049%  re¬ 
spectively.  He  reported  that  a  flexural 
strength  from  20  to  1 00  kPa  could  be  obtained, 
and  that  the  strength  measured  by  upward 
loading  was  about  70%  of  that  measured  by 
downward  loading  of  the  cantilever  beam. 
The  ratio  of  unaxial  compressive  strength 
with  horizontal  loading  to  flexural  strength 
(ojo f)  was  about  3,  while  with  vertical  load¬ 
ing  this  ratio  was  on  the  order  of  5.  The  ratio 
of  strain  modulus  E  (plate  deflection  method) 
to  Of  was  found  to  fall  within  the  range  1500- 
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2500,  similar  to  carbamide  ice  of  same  thick-  tests  in  ice  sheets  whose  thickness  and  flexu- 
ness.  Finally,  the  Kic  of  the  EG/AD/S  ice  ap-  ral  strength  are  scaled  by  the  geometric  scale 

pears  to  be  much  lower  than  that  of  carbamide  X.  The  other  properties  are  then  only  meas- 

ice  and  better  simulates  the  critical  stress  in-  ured  and  recorded.  In  all  cases  it  is  impor- 

tensity  factor  of  sea  ice  over  the  usual  range  tant  not  only  to  report  the  measured  values  of 

of  scaling  factors.  these  properties  but  also  their  measurement 

From  the  above  results,  it  appears  that  EG/  techniques  if  meaningful  comparisons 

AD/S  ice  is  a  significant  improvement  over  among  test  data  are  to  be  attempted, 

urea  ice  as  a  modeling  material.  So  far  it  has  Flexural  strength  is  usually  measured 

been  adopted  only  by  the  Institute  for  Marine  from  in  situ  tests  of  small  cantilever  beams 

Dynamics  in  St.  John’s,  Newfoundland.  Be-  and  the  effective  or  strain  modulus,  E,  by  the 

cause  of  its  cost  (about  twice  that  of  urea  ice)  plate  method.  The  kinetic  friction  coefficient 

and  the  time  involved  in  fine-tuning  the  is  measured  by  pulling  a  loaded  ice  sample 

growth  and  testing  procedures  for  each  partic-  over  a  flat  area  of  the  ship  hull  or  a  test  board 

ular  ice  tank,  other  organizations,  which  had  whose  surface  was  prepared  in  the  same 

recently  adopted  carbamide  ice,  are  awaiting  manner  as  the  ship  hull  (see  Fig.  1).  The 

or  conducting  independent  studies  of  this  new  pulling  force  is  measured  and  divided  by  the 

model  ice  before  deciding  whether  to  adopt  it  total  normal  load  to  yield  the  friction  coeffi- 

or  not.  It  can  be  mentioned  that  the  St.  John’s  cient.  Tests  may  be  conducted  at  different 

tank  initially  experienced  severe  growth  of  pulling  speeds  and  normal  pressures  to  deter- 

bacteria.  This  bacteria  problem  was  resolved  mine  possible  effects  of  these  two  parameters 

by  deleting  the  sugar  from  the  dopants  and  on  f  (ITTC  1984).  Test  methods  for  other  prop¬ 
installing  a  sand-type  filtering  system  for  erties  may  vary  greatly  from  one  test  facility 

the  tank  water.  to  another  when  those  properties  are  meas¬ 

ured  at  all.  For  refrigerated  facilities,  the 
size  of  the  tank  and  its  refrigeration  charac- 
MODEL  TEST  PROCEDURES  teristics  will  dictate  the  time  for  ice  to  grow  to 

a  specified  thickness  and  the  tempering  time 
Ice  growth  and  monitoring  required  to  achieve  a  particular  flexural 

Over  the  usual  full  range  of  geometric  strength.  When  the  target  flexural  strength  is 

scale,  X  =  10  to  40,  no  existing  model  ice  can  nearly  reached,  there  is  usually  only  a  nar- 

scale  all  the  properties  of  sea  or  freshwater  row  time  window  during  which  the  ice  proper- 

ice.  It  is  then  customary  to  conduct  ship  model  ties  must  be  measured  and  the  tests  per- 
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a.  Diagram  of  friction  test  apparatus. 

Figure  1.  Measurement  of  ice-hull  friction  coefficient. 


b.  Friction  test  apparatus  in  operation. 


Figure  1  (cont'd).  Measurement  of  ice-hull  friction  coefficient. 


formed.  This  lack  of  available  time  often 
precludes  measuring  all  ice  properties  other 
than  <xf  and  E,  which  are  considered  as  the 
most  important  to  measure,  unless  sufficient 
personnel  are  available  to  conduct  the  meas¬ 
urements  simultaneously  or  nearly  so.  To 
add  to  this  problem,  the  weaker  the  ice,  the 
narrower  the  time  window  available  and  the 
more  difficult  and  time-consuming  the 
measurements. 


EHP  tests 

Resistance  or  EHP  (effective  horsepower) 
tests  in  level  ice  are  conducted  in  a  manner 
similar  to  traditional  resistance  tests  in 
calm  water.  The  model  hull  is  connected  to  a 
towing  carriage,  free  to  heave,  pitch,  roll  and 
possibly  yaw.  Depending  on  whether  the  tow¬ 
ing  mechanism  is  stiff  or  soft,  the  hull  may 
be  totally  or  only  partially  restricted  in  surge 
and  sway.  One  major  difference  between 
open  water  and  ice  resistance  tests  is  that  in 
the  latter  the  hull  motion  is  only  steady  in  the 
mean.  As  the  hull  penetrates  into  and  crushes 
and  bends  the  ice,  the  resistance  increases 
until  the  ice  fails,  at  which  point  it  decreases 
suddenly.  This  quasi-periodic  or  cyclic  phe¬ 
nomenon,  an  example  of  which  is  shown  in 
Figure  2,  requires  that  digital  sampling  of  the 


analog  signals  be  done  at  high  enough  rates 
to  avoid  aliasing  effects.  Once  the  towing 
carriage  has  reached  a  steady  speed,  the  mod¬ 
el  resistance  is  averaged  over  a  sufficient 
travel  distance  (one  model  length  or  more)  to 
yield  reliable  values.  Since  model  tests  in  ice 
are  far  more  expensive  than  those  in  ice-free 
water,  as  many  tests  as  possible  are  to  be  con¬ 
ducted  in  one  ice  sheet. 

In  addition  to  resistance,  the  most  common 
quantities  measured  during  EHP  tests  are 
heave,  pitch  and  roll. 

The  tests  results  are  sometimes  presented 
as  the  total  resistance,  R,  versus  velocity,  V. 
However,  the  whole  range  of  tests  usually  re¬ 
quires  more  than  one  ice  sheet,  and  the  ice 
characteristics  —  thickness,  strength,  strain 
modulus,  etc.  —  will  vary  slightly  from  one 
ice  sheet  to  the  next.  This  writer  prefers  to  plot 
the  data  in  dimensionless  form  as  shown  in 
Figure  3. 

In  addition  to  quantitative  measurements, 
qualitative  test  documentation  by  means  of 
photograph,  films  or  video  recording  are  ex¬ 
tremely  useful.  In  particular,  underwater 
photography  or  video  showing  the  motion  of 
the  ice  floes  along  the  mo^e1  hull  can  be  used 
to  indicate  potential  interference  between  ice 
and  appendages.  These  visual  aids  may  sug- 
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Figure  2.  Example  of  resistance  signal  in  level  ice. 


Figure  3.  Example  of  nondimensional  resistance 
test  results. 


gest  modifications  to  hull  design  and  append¬ 
age  location  to  minimize  such  interaction. 

SHP  tests 

As  in  open  water  tests,  propulsion  or  SHP 
(shaft  horsepower)  tests  in  ice  can  be  conduct¬ 
ed  either  with  a  towed  model  or  with  a  free 
model  that  is  tethered  to  the  towing  carriage 
only  by  the  power  and  analog  data  signal 
cables.  In  the  former  method,  two  to  three  tests 
at  different  propeller  speeds  for  each  model 
speed  must  be  conducted.  The  force  (pull)  be¬ 
tween  carriage  and  model  must  be  measured. 


and  the  test  data  are  interpolated  to  zero  pull  to 
determine  the  self-propulsion  point  (rpm, 
thrust,  torque)  at  the  particular  speed.  In  the 
free-model  tests,  the  propeller  speed  is  usual¬ 
ly  set  at  a  particular  value,  and  the  corre¬ 
sponding  model  velocity  is  then  measured, 
together  with  propeller  torque  and  thrust. 
Since  model  speed  in  ice  is  only  quasi¬ 
steady,  this  method  requires  that  the  main 
carriage  be  able  to  follow  the  model  within  a 
fairly  close  distance,  and  model  speed  or  po¬ 
sition  relative  to  the  main  carriage,  or  both, 
must  be  measured.  This  test  method  has  the 
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advantage  of  yielding  the  self-propulsion 
points  without  any  interpolation  or  other  ad¬ 
justments  of  the  test  data.  On  the  other  hand, 
the  method  requires  a  rather  long  tank,  espe¬ 
cially  at  high  speed,  to  allow  the  model  to 
reach  its  quasi-steady  state,  and  is,  therefore, 
more  appropriate  to  the  longer  ice  basins. 

Since  one  of  the  goals  of  propulsion  tests  in 
ice  is  to  investigate  potential  ice-propeller 
interaction,  high-speed  data  acquisition  sys¬ 
tems  are  required  to  sample  the  thrust  and 
torque  signals  at  high  enough  rates  (about  10 
times  per  revolution). 

Examples  of  propulsion  tests  results  under 
the  captive-model  method  are  shown  in  Fig¬ 
ure  4.  As  for  resistance  tests,  visual  docu¬ 
mentation  is  very  useful,  especially  of  any 
ice-propeller  interaction. 

TEST  DATA  ANALYSIS  — 

COMPARISON  WITH  FULL  SCALE 

Analysis  of  test  results 

Ship  model  tests  are  conducted  to  predict  ei¬ 
ther  the  power  required  to  achieve  design  per¬ 
formance  or  the  performance  to  be  expected 
for  a  given  shaft  power.  In  the  first  case,  pre¬ 
dictions  are  usually  plotted  as  power  versus 
speed  in  level  ice  for  several  ice  thicknesses. 
In  the  latter  case,  the  expected  speed  at  full 
power  is  plotted  versus  ice  thickness. 

Each  ice  testing  facility  has  its  particular 
technique  for  analyzing,  interpreting  and 
extrapolating  the  test  data  to  full-scale  condi¬ 
tions.  One  technique  is  to  perform  a  regres¬ 
sion  analysis  on  the  raw  test  data.  For  exam¬ 
ple,  the  measured  ice  resistance  is  assumed  to 
be  a  polynomial  function  of  speed,  ice  thick¬ 
ness,  and  ice  flexural  strength.  Another 
method  is  to  break  the  total  resistance  into 
several  components  and  perform  regression 
analyses,  assuming  different  forms  of  the  re¬ 
gression  equation  for  each  component.  Still 
another  technique  is  to  make  similar  analy¬ 
ses  in  a  dimensionless  form. 

Extrapolation  of  the  model  test  data  to  full- 
scale  conditions  is  made  either  strictly  by  us¬ 
ing  the  Froude  scaling  law  or  by  also  intro¬ 
ducing  empirical  correction  factors  based  on 
past  experience  or  on  available  full-scale 
data  (see,  for  example,  Enkvist  and  Musta- 
maki  1986). 


In  all  cases,  because  of  the  scatter  inherent 
in  test  data  from  ice  model  tests,  seme  sort  of 
regression  analysis  is  used  to  smooth  the  re¬ 
sults.  The  analysis  may  be  made  on  the  mod¬ 
el  data  or  the  extrapolated  values,  in  dimen¬ 
sional  or  dimensionless  form.  Based  on  this 
author’s  experience,  and  until  a  more  reli¬ 
able  theoretical  analysis  of  ship-ice  interac¬ 
tion  is  available,  the  form  of  the  regression 
should  be  kept  as  simple  as  possible.  It  should 
be  selected  by  visual  inspection  of  plots  of  the 
dependent  variables  versus  the  independent 
variables  that  are  believed  to  be  of  primary 
importance.  Since  any  regression  analysis 
is  valid  only  over  the  range  of  test  parameters 
(or  the  corresponding  full-scale  values),  the 
final  results  should  not  be  extrapolated  be¬ 
yond  that  range,  which  needs  to  be  specified. 

In  the  case  of  resistance  tests,  the  data 
analysis  method  used  by  this  author  consists 
first  of  subtracting  the  open  water  resistance 
from  the  total  resistance  in  level  ice  to  obtain 
the  net  ice  resistance.  The  open  water  resis¬ 
tance  is  measured  by  towing  the  model  in  the 
ice  model  basin  free  of  ice.  The  dimension¬ 
less  net  ice  resistance  is  plotted  versus  the 
Froude  number  with  the  dimensionless  ice 
strength  as  a  parameter.  A  regression  analy¬ 
sis  is  performed  on  the  basis  of  these  plots,  the 
result  of  which  is  applied  to  full-scale  condi¬ 
tions  to  predict  the  ice  resistance  of  the  proto¬ 
type.  The  open  water  resistance  of  the  proto¬ 
type,  determined  from  traditional  open  water 
tests,  is  added  to  the  ice  resistance  to  give  the 
total  resistance  in  level  ice  at  full  scale. 
While  the  open  water  resistance  is  usually 
small  compared  to  the  ice  resistance,  in  thin 
ice  at  high  speed  (6  to  8  knots  [3  to  4  m/s])  it 
can  become  a  sizable  component  of  the  total 
resistance.  The  results  of  EHP  and  SHP 
model  tests  are  used  to  evaluate  the  thrust  de¬ 
duction  factor,  t,  and  propulsive  efficiency, 
Tjp,  as  functions  of  the  advance  coefficient  J  = 
v/nd.  Until  better  understanding  of  ice- 
propeller  interaction,  and  possible  corre¬ 
sponding  scale  effects,  is  reached,  the  model 
results  for  t{J)  and  t ;p(J)  are  assumed  to  be  ap¬ 
plicable  to  full  scale. 

Examples  of  full-scale  performance  pre¬ 
dictions  following  the  above  technique  of  data 
analysis  are  shown  in  Figure  5. 
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Figure  5.  Full-scale  resistance  prediction  from 
model  test  results. 


Comparison  with  full-scale  data 

The  validity  of  model  test  results  can  be  as¬ 
certained  only  by  comparison  with  full-scale 
trial  measurements.  There  are  very  few  ice¬ 
breakers  and  other  ice-going  vessels  that 
have  been  tested  and  whose  full-scale  trials 
results  are  not  proprietary.  Even  among  the 
full-scale  data  available  (e.g.,  Vance  1980, 
Vance  et  al.  1981,  Edwards  et  al.  1981,  Mi- 
chailidis  and  Murdey  1981)  there  are  often 
gaps  in  the  data  presented,  or  the  full-scale 
conditions  are  significantly  different  from 
the  scaled  up  model  test  values.  This  is  not  al¬ 
together  surprising  when  one  considers  the 
difficulties  involved  in  full-scale  ice  trials. 
Furthermore,  it  is  extremely  rare  to  encoun¬ 
ter  in  the  field  level  ice  sheets  devoid  of  snow, 
while  the  model  ice  sheets  are  always  so.  In 
addition,  the  presence  of  snow  may  signifi¬ 
cantly  increase  the  full-scale  resistance,  de¬ 
pending  upon  the  thickness  and  consistency 
(dry  or  wet  snow)  of  the  snow  cover. 

In  this  author’s  opinion,  some  of  the  best 
sets  of  full-scale  data  currently  available  are 
those  collected  by  Vance  (1980)  and  Vance  et 
al.  (1981)  with  the  USCGC  Katmai  Bay,  and 
by  Michailidis  and  Murdey  (1981)  with  the 
CCGC  Franklin.  Comparison  of  the  ship  per¬ 
formance  predictions  from  model  tests  con¬ 
ducted  at  CRREL  (Tatinclaux  1984a.b,  1985) 


with  full-scale  measurements  for  both  of 
these  ships  are  shown  in  Figure  6.  As  can  be 
seen,  the  power  predictions  from  model  tests 
are  consistently  greater  than  field  measure¬ 
ments.  This  author  believes  that  this  discre¬ 
pancy  is  at  least  in  part  attributable  to  incor¬ 
rect  modeling  of  the  ice  piece  size  (Keinonen 
1983)  resulting  from  incorrect  strain  modu¬ 
lus,  fracture  toughness  and  density  of  the 
model  ice.  This  may  lead  to  exaggerated,  un¬ 
realistic  ice-propeller  interaction  in  the 
model. 


ANALYTICAL  AND  EMPIRICAL 
PREDICTORS 

In  parallel  with  the  development  of  labora¬ 
tory  model  testing  techniques,  efforts  have 
been  made  to  develop  predictor  equations  or 
numerical  schemes  for  the  level  ice  resis¬ 
tance  and  propulsion  performance  of  ice¬ 
going  vessels.  These  equations  or  schemes 
can  be  divided  into  two  categories:  1)  empiri¬ 
cal  and  2)  analytical  or  semi-analytical. 

Empirical  predictors 

Empirical  predictor  equations  for  the  re¬ 
sistance  or  propulsion  characteristics  in  lev¬ 
el  ice,  or  both,  have  been  proposed  by  Kashtel- 
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a.  USCGC  Katmai  Bay. 
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Figure  6.  Comparison  between  predictions  and 
full-scale  trials  measurements. 


yan  et  al.  (1968),  Enkvist  (1972),  Lewis  et  ai. 
(1982),  and  Tsoy  (1983).  The  equations  were 
obtained  by  multiple  regression  analysis  of 
available  full-scale  and  model  test  data.  As 
for  any  empirical  formula,  these  equations 
are  valid  only  over  the  range  of  parameters 
for  which  they  were  determined.  In  particu¬ 
lar,  the  data  bases  for  all  the  equations  were 
obtained  from  classical  icebreaker  hull 
shapes.  Therefore,  it  is  unlikely  that  the  pre¬ 
dictors  may  be  used  for  the  new  hull  shapes 
currently  under  development  and  discussed 
in  a  following  section. 

The  set  of  equations  proposed  by  Tsoy 
(1983)  was  found  by  this  writer  to  be  of  particu¬ 
lar  practical  interest.  It  predicts  the  limiting 
ice  thickness  capable  of  being  broken  by  a 
given  ship,  i.e.,  her  design  point  in  level  ice. 
Given  the  total  shaft  power,  Pw  in  kilowatts, 
and  the  propeller  diameter,  efp  in  meters,  the 
thrust  (Th)  developed  at  near-mooring  operat¬ 
ing  conditions  is  estimated  by 

Th  =  *i  (rfp  •  Pw)  2/3  (kN) 

where  is  an  empirical  coefficient:  k1  =  0.78 
for  a  single-screw  ship,  kY  =  0.98  for  a  double¬ 
screw  ship,  and  k j  =  1.12  for  a  triple-screw 
ship,  assuming  all  propellers  to  be  identical 
and  that  power  is  equally  distributed  among 
the  propellers. 

The  limiting  ice  thickness,  hi  in  meters, 
that  the  ship  is  capable  of  breaking  in  the  con¬ 
tinuous  mode  (i.e.,  at  a  speed  of  3  knots  [1.5 
m/s|)  is  then  given  by 


h^  —  k  2 A 

where  B  is  the  ship  beam  in  meters,  A  is  the 
ship  displacement  in  metric  tonnes,  and  k2  is 
an  empirical  coefficient  given  by  Tsoy  as 

k2  =  0.031  to  0.035 

for  icebreakers,  and 

k2=  0.039  -0.025 P 

for  ice-reinforced  ships,  with  /?  being  the  ship 
block  coefficient. 


Once  this  limiting  ice  thickness  is  esti¬ 
mated,  and  when  the  maximum  ship  speed  in 
open  water  is  known,  the  ship  performance 
curve,  maximum  speed  vs  ice  thickness,  can 
be  drawn  under  the  assumption  that  it  is  line 
ar.  Such  performance  curve  was  estimated 
for  the  USCGC  Katmai  Bay  and  compared 
with  model  tests  predictions  in  Figure  7. 

Analytical  and  semi-analytical  schemes 

All  analytical  and  semi-analytical  deri 
vations  of  ship  resistance  in  level  ice  (Mila 
no  1973,  Carter  1982,  Kotras  et  al.  1983;  avail 
able  in  the  literature  assume  that  the  ice 
resistance  can  be  divided  into  several  compo¬ 
nents,  namely  a  breaking  resistance  due  to 
ice  crushing,  shearing,  bending  or  buckling 
(or  all  four),  a  frictional  resistance,  and  a 
component  due  to  ice  floe  submergence  and 
entrainment.  The  various  authors  follow  ei¬ 
ther  a  force  balance  (Carter,  Kotras  et  al.)  oi 
an  energy  balance  (Milano)  approach.  The 
mathematical  models  of  Milano  and  Carter 
appear  to  be  purely  analytical,  since  the  only 
inputs  required  are  the  ice  properties  and  the 
ship  hull  geometry.  The  resistance  equation 
proposed  by  Kotras  et  al.  is  considered  semi- 
analytical,  since  it  still  contains  empirical 
coefficients  that  need  to  be  adjusted  for  each 
ship  studied  from  a  few  model  test  data  or 
full-scale  trial  data;  the  equation  may  then  be 
applied  to  conditions  other  than  those  used  in 
determining  the  coefficients. 

Milano  has  applied  his  method  to  numer¬ 
ous  ships  of  various  shapes  and  sizes  (e.g., 
Milano  1982)  with  reasonable  correlation  be¬ 
tween  predicted  and  measured  ice  resistance. 
Milano  also  claims  that  his  analytical  model 
is  applicable  to  any  hull  form.  Carter’s  (1982) 
comparisons  between  measured  and  predict¬ 
ed  resistance  are  also  quite  satisfactory,  in 
spite  of  his  equation  totally  neglecting  iner¬ 
tial  and  buoyancy  forces  due  to  submergence 
of  ice  floes  by  the  moving  vessel.  Further¬ 
more,  his  treatment  of  the  effect  of  ship  veloci¬ 
ty  on  ice  resistance  remains  unconvincing  to 
this  writer.  Finally,  his  equation  is  based  in 
part  on  elliptical  approximations  of  the  ship 
waterlines,  and  therefore  may  not  be  readily 
applicable  to  the  new,  unconventional  ice¬ 
breaker  hull  forms  currently  under  develop¬ 
ment. 
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Figure  7  USCGC  Katmai  Bay  performance  curve  in 
level  ice  measured  during  trials,  predicted  from  model 
tests,  calculated  by  Tsoy’s  (1984)  method. 


CURRENT  RESEARCH  EFFORTS 
IN  ICE  MODELING 

Research  and  development  efforts  at  ice 
modeling  facilities  are  pursued  in  many  di¬ 
rections:  improvement  of  model  ice,  field 
and  laboratory  measurements  or  characteri¬ 
zations  of  ice  properties,  scale  effects  and  test 
data  corrections,  ice-propeller  interaction, 
etc.  International  cooperative  efforts  between 
research  facilities  are  also  underway  under 
the  aegis  of  the  ITTC  Committee  on  Perfor¬ 
mance  of  Ships  in  Ice-Covered  Waters  (Ice 
Committee  for  short). 

At  the  outset  it  can  be  pointed  out  that  be¬ 
cause  of  particular  characteristics  and  condi¬ 
tions  (size,  refrigeration  system  type  and  ca¬ 
pacity,  local  environment,  etc.)  each  ice  tank 
has  developed  its  own  method  for  manufac¬ 
turing  a  level  ice  sheet  of  prescribed  mechan¬ 
ical  properties.  Modifications  to  seeding, 
growth  monitoring  and  tempering  tech¬ 
niques  (in  the  case  of  doped  ice)  are  constant¬ 
ly  introduced  to  improve  the  model  ice  consis¬ 
tency  and  repeatability. 


International  cooperative  research 

One  of  the  tasks  of  the  ITTC  Ice  Committee 
is  to  develop  standard  methods  for  measure¬ 
ment  of  ice  properties  both  in  the  field  and  in 
the  laboratory.  It  is  also  attempting  to  develop 
guidelines  showing  which  quantities  need  to 
be  documented  during  model  tests  and  full- 
scale  trials.  The  purpose  of  these  standards 
and  guidelines  is  to  allow  meaningful  com¬ 
parison  between  results  obtained  at  different 
facilities.  Obviously,  they  need  to  be  periodi¬ 
cally  revised  as  progress  in  the  understand¬ 
ing  of  ship-ice  interaction,  testing  methodol¬ 
ogy,  and  testing  equipment  is  made. 

The  main  recent  effort  of  the  Ice  Committee 
was  the  testing  of  a  standard  model  by  all 
participating  ice  tanks.  The  standard  model 
selected  was  that  of  the  Canadian  Coast 
Guard  “R”-class  icebreaker.  Fiberglas  mod¬ 
els  at  the  scales  of  1:20  and  1:40  were  con¬ 
structed  by  the  National  Research  Council  of 
Canada  and  circulated  to  the  participants. 
The  ice  research  facilities  who  participated 
in  the  joint  research  program  were  IMD/NRC 
(Canada),  CRREL  (USA),  HSVA  (Germany), 


19 


WARC  (Finland),  AARI  (USSR),  JSRI  and 
NKK  (Japan).  The  test  program  (ITTC  1981) 
included  resistance  and  propulsion  tests  at 
two  ice  thicknesses  and  two  ice  strengths. 
The  results  of  the  first  series  of  tests  were  re¬ 
ported  to  the  International  Towing  Tank 
Conference  in  Goteborg,  Sweden  (ITTC  1984). 
This  international  cooperative  research  pro¬ 
gram  is  continuing,  the  results  of  which  are 
to  be  reported  periodically  to  the  ITTC. 

Ice  testing 

Testing  for  the  mechanical  properties  of 
ice  is  a  difficult  and  time-consuming  task. 
Furthermore,  all  current  methods  give  only 
local  values  with  the  assumption  that  these 
values  are  representative  of  the  whole  ice 
sheet.  Baker  et  al.  (1987)  have  proposed  an 
ice-indexer  for  laboratory  use  that  would 
measure  ice  flexural  strength,  shear  strength 
and  crushing  strength  simultaneously  and 
continuously.  This  ice-indexer  is  still  under 
development  and  evaluation. 

While  measurement  of  the  friction  coeffi¬ 
cient  of  ice  is  standard  procedure  in  the  labor¬ 
atory,  no  satisfactory  field  apparatus  had 
been  devised.  Kitagawa  and  Izumiyama 
(1986)  have  recently  reported  on  a  portable 
field  device  for  ice  friction  measurements 
that  could  be  a  significant  addition  to  the  tools 
needed  to  characterize  field  ice. 

It  is  only  recently  that  the  fracture  tough¬ 
ness  of  ice  has  received  increasing  attention 
as  a  material  property  significantly  affect¬ 
ing  ship  resistance  in  ice.  As  far  as  this  writ¬ 
er  knows,  the  Institute  for  Marine  Dynamics 
of  the  National  Research  Council  of  Canada, 
St.  John’s,  Newfoundland,  is  the  only  ice 
modeling  facility  where  the  fracture  tough¬ 
ness  of  its  model  ice  or  some  related  index  is 
measured  routinely.  No  publication  is  yet 
available  on  the  relation  between  fracture 
toughness  and  ship  model  resistance,  or  on 
how  to  modify  the  full-scale  prediction  to  in¬ 
clude  Kjc.  CRREL  has  recently  initiated  an 
investigation  of  the  fracture  toughness  of 
urea-ice  (Bentley  et  al.,  in  press),  a  study  to  be 
expanded  to  freshwater  ice.  Standard  meas¬ 
urement  procedures  and  apparatus  are  yet  to 
be  devised. 

CRREL  research  on  ship-ice  interaction 

As  was  mentioned  earlier,  in  this  author’s 
experience  the  power  requirements  predicted 


by  direct  Froude  scaling  of  model  tests  results 
were  consistently  greater  than  those  meas¬ 
ured  during  full-scale  trials.  From  visual 
observation  via  underwater  video  and  con¬ 
versation  with  icebreaker  operators,  in  par¬ 
ticular  personnel  of  the  U.S.  Coast  Guard,  the 
conclusion  was  reached  that  ice  ingestion  by 
the  propellers  is  far  greater  in  the  model  tests 
than  at  full  scale.  This  conclusion  is  in 
agreement  with  observations  by  Keinonen 
(1983)  who  attributed  this  excessive  ice- 
propeller  interaction  to  model  ice-floes  being 
relatively  larger  than  full-scale  floes. 

This  author  investigated  the  size  distribu¬ 
tion  of  ice  floes  created  by  a  simple  wedge 
(idealized  icebreaker  bow)  towed  through  ice 
sheets  of  varying  thickness,  flexural 
strength  and  strain  modulus  (Tatinclaux 

1986) .  It  was  found  that  the  plan  area  and 
largest  dimension  of  the  floes  followed  a  log¬ 
normal  probability  density  function  (Fig.  8). 
Similar  results  were  presented  by  McKindra 
and  Lutton  (1981)  from  trial  data  with  a  USCG 
Great  Lakes  icebreaker  (WTGB  class).  The 
ratio  of  model  floes’  mean  plan  area  to  ice 
thickness  squared,  PJh"\,  was  found  to  be  only 
proportional  to  the  dimensionless  strength, 
Cf  ly^,  but  independent  of  other  parameters 
such  as  the  strain  modulus  or  characteristic 
length  (Fig.  9).  The  coefficient  of  proportion¬ 
ality  depended  on  the  type  of  model  ice  (urea- 
doped  or  synthetic  ice).  The  available  full- 
scale  data  were  insufficient  to  confirm  these 
last  results. 

A  subsequent  study  investigated  the  effect 
of  floe  size  on  propeller  torque  (Tatinclaux 

1987) .  Captive  model  propulsion  tests  were 
made  with  the  1:10  model  of  the  USCGC  Kat- 
mai  Bay  previously  tested  at  CRREL.  The 
tests  were  made  in  pre-sawn  channels  with 
floe  sizes  ranging  from  1/2  to  1/6  the  model 
beam,  and  in  brash-ice-filled  channels.  At  a 
given  advance  coefficient,  the  torque  coeffi¬ 
cient  E"q,  increased  linearly  with  floe  size 
(Fig.  10),  and  the  torque  standard  deviation 
Sq  increased  with  velocity,  the  rate  of  in¬ 
crease  being  some  function  of  floe  size.  Also, 
it  appeared  that  below  a  certain  ship  speed  no 
ice-propeller  interaction  occurred,  and  that 
this  critical  velocity  decreased  with  increas¬ 
ing  floe  size  (Fig.  11). 

On  the  basis  of  the  above  study,  it  is  suspect¬ 
ed  that  model  ice  specific  gravity  and  the  ice- 
hull  friction  coefficient  may  significantly 
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Figure  1 1.  Effect  of  floe  size  and  velocity 
on  propeller  torque  standard  deviation 
(Tatinclaux  1987). 


iffect  the  trajectory  of  the  ice  floes  along  the 
mil  and  therefore  the  ice-propeller  interac- 
ion.  As  a  ship  (model  or  prototype)  travels 
h rough  ice,  the  broken  floes  are  submerged 
)y  the  bow  to  some  depth,  and  they  then  rise  be¬ 
cause  of  buoyancy.  Ice  may  be  ingested  if  the 
;hip  has  traveled  one  ship-length  before  the 
loes  have  fully  risen  to  near  the  surface. 
Since  model  ice  is  usually  somewhat  denser 
.han  field  ice,  the  model  floes  may  be  sub- 
nerged  deeper  and  rise  slower  than  full- 
icale  ice  floes.  Therefore,  it  can  be  expected 


that  ice  will  be  ingested  by  the  propeller  at  a 
ship  speed  lower  than  at  full  scale,  and  that  it 
will  also  be  more  severe.  An  analytical  and 
experimental  study  has  been  initiated  to  ad¬ 
dress  these  questions.  The  results  of  a  first 
crude  mathematical  model  of  ice  floe  trajec¬ 
tory  along  a  rigid  surface  are  shown  in  Fig¬ 
ure  12. 

As  a  final  comment,  it  is  again  empha¬ 
sized  that  one  main  difficulty  in  comparing 
results  of  model  studies  with  full-scale  re¬ 
sults  is  the  lack  of  reliable  and  complete  sets 
of  prototype  data. 


NOVEL  ICEBREAKING  BOW  DESIGNS 

The  past  decade  has  witnessed  the  develop¬ 
ment  of  new  icebreaking  bows,  the  shapes  of 
which  depart  radically  from  the  conventional 
ones.  A  bow  shape,  first  proposed  and  patented 
by  Waas  (1976),  has  been  further  developed 
and  extensively  tested  in  West  Germany  both 
in  the  laboratory  and  in  the  field  (Freitas 
1981, 1982;  Hellmann  1982;  Freitas  and  Nish- 
izaki  1985;  Schwarz  1986).  From  its  original 
conception  (Fig.  13)  the  Waas  bow  has 
evolved  into  radically  new  shapes  (Fig.  14). 
Basically,  it  has  square  shoulders,  a  very  low 
stem  angle  at  the  water  line,  and  an  initially 
flat  bottom  that  evolves  into  a  wedge.  This 
bow  first  shears  the  level  ice  at  the  shoulders, 
then  breaks  it  in  bending  into  pairs  of  regu¬ 
larly  shaped  ice  floes,  which  are  deflected 
outwards  under  the  remaining  ice  sheet.  Ac¬ 
cording  to  its  developers,  the  main  advantag¬ 
es  of  the  Waas  bow  over  conventional  ice¬ 
breaker  bows  are:  more  efficient  icebreak- 


Figure  12.  Expected  effect  of  ice  density  on  ice  floe  trajectory  along 
a  ship  hull. 
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Figure  14.  Recent  developments  in  Thyssen-Waas  bow  (From 
Transactions  of  the  ASME,  Journal  of  Energy  Resources, 

A.  Freitas  and  R.S.  Nishizaki,  June  1986). 


Figure  15.  Wartsila  experimental  bow  (From  Enkvist  and  Mustamaki 
1986). 


ing,  i.e.,  reduced  ice  resistance  and  required 
power;  improved  maneuverability  in  level 
ice;  and  creation  of  a  nearly  ice-free  channel 
with  little  or  no  ice-propeller  interaction  up  to 
speeds  of  5  knots  [2.5  m/s).  The  effects  of  the 
Waas  bow  on  a  vessel’s  seakeeping  capabili¬ 
ties,  performance  in  waves  and  ice  ridges  re¬ 
main  a  subject  of  controversy.  A  soviet  ice¬ 
breaker  has  recently  been  retrofitted  with 
such  a  bow,  the  U.S.  Coast  Guard  has  initiat¬ 
ed  an  evaluation  program  of  its  performance 
potential,  and  its  developers  are  attempting  to 


adapt  this  bow  to  shallow  draft  icebreakers  for 
operation  in  estuaries  and  rivers. 

Another  experimental  icebreaking  bow 
has  been  recently  designed  and  tested  at  mod¬ 
el  and  full  scales  in  Finland  (Enkvist  and 
Mustamaki  1986).  This  bow  has  a  small  stem 
angle,  circular  waterlines  and  a  plow  or  de¬ 
flector  vanes  to  deflect  the  broken  ice  floes 
outwards  (Fig.  15).  Model  and  full-scale  tests 
showed  a  significant  decrease  in  icebreak¬ 
ing  resistance  in  level  ice  over  a  convention¬ 
al  bow,  but  no  essential  improvement  of  per- 


formance  in  ridges  and  old  broken  channels, 
and  no  improvement  in  maneuverability  in 
ice.  An  increase  in  open  water  resistance 
and  in  slamming  during  transit  in  head  sea 
was  observed.  Enkvist  and  Mustamaki  did 
not  recommend  this  innovative  bow  for  ves¬ 
sels  operating  primarily  in  open  water  or  in 
ice  ridges  or  rubble  ice.  They  would  limit  the 
operation  of  this  bow  to  areas  of  thick  level  ice 
such  as  rivers,  lakes  and  sheltered  regions 
along  a  sea  coast.  In  the  oral  presentation  of 
their  paper,  Enkvist  and  Mustamaki  indicat¬ 
ed  that  the  experimental  bow  was  mounted  on 
a  towboat  and  used  to  open  the  Saima  Canal  to 
navigation  at  the  end  of  the  1985-86  winter. 
They  reported  that  a  navigation  channel, 
nearly  free  of  ice,  was  opened  in  less  than  a 
day,  when  in  previous  years  it  had  taken  sev¬ 
eral  days  for  a  towboat  working  alone  to  open 
the  navigation  channel,  which  contained  a 
large  amount  of  brash  ice. 

While  not  radically  different  from  con¬ 
ventional  icebreaker  bows,  a  new  bow  was  re¬ 
cently  designed  by  Melville  Shipping  Ltd,  Ot¬ 
tawa,  for  the  cargo  ship  MV  Arctic.  This 
Melville  bow  proved  to  significantly  improve 
the  ice  performance  of  the  MV  Arctic  (Baker 
and  Nishizaki  1986),  confirming  the  model 
tests  carried  out  at  various  facilities. 


CONCLUSIONS  AND  FINAL  REMARKS 

Since  the  mid-1970s  significant  improve¬ 
ments  have  been  made  in  the  methodology  of 
ship  testing  in  ice  for  both  EHP  and  SHP  tests, 
as  well  as  in  the  development  and  characteri¬ 
zation  of  model  ice.  The  art  and  science  of 
model  testing  in  ice  has  now  reached  the  point 
where  model  tank  operators  are  confident 
that  model  test  results  can  be  used  in  predict¬ 
ing  the  performance  of  ice-going  vessels,  in 
suggesting  modifications  to  existing  ships  to 
improve  their  ice  performance  (e.g.,  Baker 
and  Nishizaki  1986),  and  in  designing  radi¬ 
cally  new  bow  shapes  for  icebreakers  and  oth¬ 
er  ice-going  vessels.  The  progress  thus  made 
has  allowed  the  construction  of  ice-going  ves¬ 
sels  with  improved  ice  performance  at  lower 
power  requirements.  Further  improvements 
are  to  be  expected  in  the  future  as  better  under¬ 
standing  of  the  physics  of  ship-ice  interac¬ 
tion  is  gained  and,  especially,  as  more  and 
better  field  data  are  acquired  from  full-scale 


trials.  Areas  where  improvements  are  more 
urgently  needed  are: 

1.  Knowledge  of  the  effect  of  fracture  tough¬ 
ness  on  icebreaking  resistance,  breaking 
pattern  and  ice  floe  size  distribution. 

2.  Knowledge  of  the  effect  of  ice  floe  size 
and  model  ice  density  on  floe  trajectory,  ice- 
propeller  interaction  and  propeller  perfor¬ 
mance  (thrust,  torque  and  shaft  power). 

3.  Measurement  of  ice-hull  friction  coeffi¬ 
cient  at  model  and  full  scales  and  its  effect  on 
icebreaking  resistance  and  propeller  perfor¬ 
mance. 

4.  Analytical  description  of  ship-ice  inter¬ 
action. 

5.  Detailed  documentation  of  full-scale 
trials  in  ice,  especially  of  the  ice  mechanical 
properties. 

6.  Development  of  apparatus  and  stan¬ 
dardized  methods  for  measuring  ice  proper¬ 
ties,  both  at  model  and  full  scales;  in  particu¬ 
lar,  an  apparatus  for  field  measurement  of 
ice-hull  friction  coefficients. 

While  this  paper  dealt  only  with  modeling 
of  ships  in  level  ice,  it  should  be  mentioned 
that  research  efforts  are  also  actively  pursued 
in  areas  of  ship  motion  in  brash  and  mush 
ice,  and  of  ship-ice-ridge  interaction  (ram¬ 
ming  tests)  (ITTC  1984).  The  latter  is  of  par¬ 
ticular  importance  since  extensive  fields  of 
level  ice  are  the  exception  in  the  arctic  and 
subarctic  regions,  while  ridges  are  ever 
present  and  can  seldom  be  entirely  avoided. 
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